Presnel coefficients as hyperbolic rotations 
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We describe the action of a plane interface between two semi-infinite media in terms of a transfer 
matrix. We find a remarkably simple factorization of this matrix, which enables us to express the 
Fresnel coefficients as a hyperbolic rotation. 
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I. INTRODUCTION 

Reflection and transmission at a discontinuity are 
perhaps the first wavelike phenomena that one en- 
counters in any undergraduate physics course. The 
physics underlying this behavior is well understood: mis- 
matched impedances generate the reflected and trans- 
mitted waves, while the application of the proper bound- 
ary conditions at the discontinuity provide their corre- 
sponding amplitude coefficients jij. Moreover, this gen- 
eral framework facilitates a unified treatment for all the 
types of waves appearing in Nature. 

For light waves the impedance is proportional to the 
refractive index. Accordingly, the behavior of light at 
the plane interface between two semi-infinite media are 
derived in most optics textbooks |2|, 0, . The resulting 
amplitude coefficients are described by the famous Fres- 
nel formulas. It seems almost impossible to say anything 
new about these Fresnel formulas. However, a quick look 
at the indexes of this Journal ||, among others |], 
immediately reveals a steady flow of papers devoted to 
subtle aspects of this problem, which shows that the topic 
is far richer than one might naively expect. 

In this paper we reeleborate once again on this theme. 
We present the action of any interface in terms of a trans- 
fer matrix, and we find a hitherto unsuspectedly simple 
factorization of this matrix. After renormalizing the field 
amplitudes, such a factorization leads us to introduce a 
new parameter in terms of which Fresnel formulas appear 
as a hyperbolic rotation. 

As our teaching experience demonstrates, the students 
have troubles in memorizing Fresnel formulas due to their 
fairly complicated appearance. In this respect, there are 
at least two reasons that, in our opinion, bear out the 
interest of our contribution: first, Fresnel formulas ap- 
pears in the new variables as a hyperbolic rotation that 
introduces remarkable simplicity and symmetry. Second, 
this formalism is directly linked to other fields of physics, 
mainly to special relativity, which is more than a curios- 



II. THE INTERFACE TRANSFER MATRIX 

Let two homogeneous isotropic semi-infinite media, de- 
scribed by complex refractive indices Nq and N\ , be sep- 
arated by a plane boundary. The Z axis is chosen per- 



pendicular to the boundary and directed as in Fig. 1. 

We assume an incident monochromatic, linearly polar- 
ized plane wave from medium 0, which makes an angle 
9q with the Z axis and has amplitude Eq . The electric 
field is either in the plane of incidence (denoted by super- 
script ||) or perpendicular to the plane of incidence (su- 

(r) 

perscript _L). This wave splits into a reflected wave E Q 
in medium 0, and a transmitted wave Eq in medium 1 
that makes and angle d\ with the Z axis. The angles of 
incidence 0q and refraction B\ are related by the Sncll's 
law 



TVosinflo = iVisinfli. 



(1) 



If media and 1 are transparent (so that No and Ni are 
real numbers) and no total reflection occurs, the angles 9q 
and 0\ are also real and the above picture of how a plane 
wave is reflected and refracted at the interface is simple. 
However, when either one or both media is absorbing, 
the angles 6q and 9\ become, in general, complex and 
the discussion continues to hold only formally, but the 
physical picture of the fields becomes complicated jllj . 

We consider as well another plane wave of the same 
frequency and polarization, and amplitude e[ ' , incident 
from medium 1 at an angle 9\, as indicated in Fig. 1. In 

(r) (t) 

the same way, we shall denote by E\ and E\ the re- 
flected and transmitted amplitudes of the corresponding 
waves. 

The complex amplitudes of the total output fields at 
opposite points immediately above and below the inter- 
face will be called Eq ' and , respectively. The wave 
vectors of all waves lie in the plane of incidence and when 
the incident fields are || or _L polarized, all plane waves 
excited by the incident ones have the same polarization. 

The amplitudes Eq ^ and E± are then given by 
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where vqi and £oi ar e the Fresnel reflection and trans- 
mission coefficients for the interface 01, and rio and t\o 
refer to the corresponding coefficients for the interface 10. 
These Fresnel coefficients are determined by demanding 
that across the boundary the tangential components of E 
and H should be continuous B . For nonmagnetic media 
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they are given by 



b ax ~ 



' 01 



'01 



Ni cos 6 — N cos 9\ 
Nx cos 9 + N cos 0x 1 

2No cos O 
N x cos O + N cos 6»i ' 



iVo cos 6>o — iVi cos 6*i 
iV cos O + N x cos 6»i ' 

2iV cos O 
N cos 6» + Nx cos 0i ' 



(3) 



(4) 



represents the matrix of a clockwise rotation of angle 7r/4. 
Now, it is straightforward to convince oneself that a diag- 
onal matrix postmultiplied by R(ir/4) and premultiplied 
by its inverse is always of the form 



a b 
b a 



(9) 



a and 6 being, in general, complex numbers. This result 
implies that loi in Eq. (^) must be also of this form, which 
turns out the constraints 



no = -rax, 



1 + nnrio = ioi*io- 



(10) 



for both basic polarizations. It is worth noting that, 
although these equations are written for electromag- 
netic waves, it is possible to translate all the results for 
particle-wave scattering, since there is a one-to-one cor- 
respondence between the propagation in an interface be- 
tween two media of electromagnetic waves and of the non- 
relativistic particle waves satisfying Schrodinger equa- 
tion @. 

The linearity revealed by Eqs. (U) suggests the use of 
2x2 matrix methods. However, Eqs. (g) links output 
to input fields, while the standard way of treating this 
topic is by relating the field amplitudes at each side of 
the interface. Such a relation is expressed as |l^, [l4| 



(5) 



The choice of these column vectors is motivated from 
the optics of layered media, since it is the only way of 
calculating the field amplitudes at each side of every layer 
by an ordered product of matrices. 

We shall call loi the interface transfer matrix and, from 
Eqs. (0), is given by 




This applies to both basic polarizations by the simple at- 
tachment of a label to all the coefficients and constitutes 
an alternative algebraic demonstration of the well-known 
Stokes relations without resorting to the usual time- 
reversal argument [lq] . Similar results can be also 
derived in particle scattering from the unitarity require- 
ment on the S matrix. However, note that the equality 
\ r ia\ = l r oi| implied by Eq. (|10| ) can become counter- 
intuitive when applied to particle reflection, since one 
might expect stronger reflection for particle waves mov- 
ing up in a potential gradient than for those going down. 
In fact, these relations, as emphasized by Lekner |jl2| , 
ensure that the reflectivity is exactly the same in the two 
cases, unless there is total internal reflection. 

In summary, these Stokes relations allows one to 
write 



, W 1 r 01 

loi -%Ui 1 



(11) 



It is worth noting that the inverse matrix satisfies 
Iq] 1 = I io and then describes the interface taken in the 
reverse order. The physical meaning of these matrix ma- 
nipulations is analyzed in Section M. 
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By using a matrix formulation of the boundary condi- 
tions |l5| ] one can factorize the interface transfer matrix 
loi in the new and remarkable form (that otherwise 
one can also check directly using the Fresnel formulas) 



. ,, [ , / cos Ox I cos#n , 

'oi = R ( 17 Ni/Nq ) R(tt/4), 



(7) 



,± = R-i(7r/4)( (^icos0i)/(iVocos0o) I) 



JW/4), 



where 



W4) = 7l(i1 



(8) 



III. RENORMALIZATION OF FIELD 
AMPLITUDES 

From Eqs. (0) one directly obtain that, for both basic 
polarizations, we have 



, .11 nil -^1 COS 01 , 

de <= det '°-A^^ L 



(12) 



For the reasons that will become clear in Section [V, it 



is adequate to renormalize the field amplitudes to ensure 
that the transfer matrix has always unit determinant. To 
this end, let us define 



= V^ocosflo 4 , 
e[ ±] = y/NxCOS0x E[ ±] . 



(13) 
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Accordingly, the action of the interface is described now 
by 



J+) 
-o 
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\ °0 / \ °1 

where the renormalized interface matrix is 



1 01 



= I ^oi + 1/Coi Coi - 1/Coi 
2 " Vfri £01 + VCoi 

and the factor £oi has the values 



(14) 



(15) 



The Fresnel coefficients can be obtained from the 
renormalized ones as 



roi 



roi, 



(21) 



N x cos 9 1 
N cos 0o 



to 



It is clear from Eqs. (|l^) that the single parameter 
£oi gives all the information about the interface, even for 
absorbing media or when total reflection occurs. We have 
'oi = 'loi that is, the inverse also describes the interface 
taken in the reverse order. Thus, £io = l/£oi arLC1 h 
follows that 



*01 



iViCos6» _ /sin(20 o ) 
iY cos6>i ~~ V sin(26»i)' 



fife = 



N cos 9 I tan 9\ 
N\ cos 9\ V tan 9q 

Other way of expressing these relations is 



ell t± _ 
?01?01 — 



Coi 



COS Vo 

cos 9\ ' 
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No 



(16) 



(17) 



It is now evident from Eq. ( |15| ) that the renormalized 
interface matrix satisfies det i i = +1, as desired. More- 
over, by taking into account the general form given in 
Eq. (p"l|), we can reinterpret ioi in terms of renormalized 
Fresnel coefficients as 



where 



1 01 



1 ( 1 fbi 



i V r oi 1 J' 



goi ~ 1/goi 
&i + l/Coi 



Coi + 1/Coi 



(18) 



(19) 



which satisfy 



r Q i 2 +to! 2 = 1. 



(20) 



This relation does not trivially reduce to the conservation 
of the energy flux on the interface, because the complex 
reflection and transmission coefficients appear in the form 
f i 2 and t i 2 instead of |r i| 2 and \toi\ 2 - In fact, it can 
be seen as a consequence of the renormalization factors 
appearing in the definition ( |l3| ) that project the direction 
of the corresponding wave vector onto the normal to the 
boundary. 



roi 



toi — t 



-no, 



(22) 



lo- 



in Fig. 2 we have plotted the behavior of £q X and £qi as 
a function of the angle of incidence 9o , for an interface air- 
glass (No/Ni = 2/3) and, for the purpose of comparison, 
the corresponding values of r 01 and r^. The discussion 
about these amplitude coefficients and the corresponding 
phase shifts can be developed much in the same way as 
it is done in most of the undergraduate optics textbooks. 



IV. THE INTERFACE AS A HYPERBOLIC 
ROTATION 

The definition of the renormalized transfer matrix for 



an interface in Eq. (15) may appear, at first sight, rather 
artificial. In this Section we shall interpret its meaning 
by recasting it in an appropriate form that will reveal the 
origin of the rotation matrices R(ir/4). 

To simplify as much as possible the discussion, let us 
assume that we are dealing with an interface between 
two transparent media when no total reflection occurs. In 
this relevant case, the Fresnel reflection and transmission 
coefficients, and therefore £oi; are real numbers. Let us 
introduce a new parameter £oi by 

£oi = exp(Coi/2). (23) 

Then, the action of the interface can be expressed as 



'0 

o(-) 

"0 



cosh(Coi/2) sinh(Coi/2) 
sinh(Coi/2) cosh(Coi/2) 




(24) 



where the renormalized Fresnel coefficients can be writ- 
ten now as 



foi = tanh(Coi/2), 



(25) 
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cosh(Coi/2)' 
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Given the importance of this new reformulation of the 
action of an interface, some comments seem pertinent: 
it is clear that the reflection coefficient can be always 
expressed as a hyperbolic tangent, whose addition law is 
simple. In fact, such an important result was first derived 
by Khashan |U| and is the origin of several approaches 
for treating the reflection coefficient of layered struc- 
tures [|l7j, including bilinear or quotient functions p8[ , 
that are just of the form (|T9j). However, the transmis- 
sion coefficient for these structures seems to be (almost) 
safely ignored in the literature, because it behaves as a 
hyperbolic secant, whose addition law is more involved. 

Now the meaning of the rotation R(7r/4) can be put 
forward in a clear way. To this end, note that the trans- 
formation (24) is formally a hyperbolic rotation of angle 
Coi/2 acting on the complex field variables [e^ + \e^ - ^]. 
As it is usual in hyperbolic geometry |jl9| , it is conve- 
nient to study this transformation in a coordinate frame 
whose new axes are the bisecting lines of the original one. 
In other words, in this frame whose axes are rotated 7r/4 
respect to the original one, the new coordinates are 



g(+) 

gc-) 



= R(tt/4) 



;(+) 
J-) 



(26) 



for both and 1 media, and the action of the interface is 
represented by the matrix 



(27) 



which is a squeezing matrix that scales down to the 
factor £oi and e[ up by the same factor. 

Furthermore, the product of these complex coordinates 
remains constant 




s(+)s(-) _ s(+)s(-) 



n 



or 



,(+)l2 



'0 



,(")l2 



,(+)l2 



,(")l2 



(28) 



(29) 



which appears as a fundamental invariant of any inter- 
face. In these renormalized field variables it is nothing 
but the hyperbolic invariant of the transformation. When 
viewed in the original field amplitudes it reads as 

NoCoseoUE^] 2 -^] 2 } = iVrCOS^W] 2 -^] 2 }, 

(30) 

which was assumed as a basic axiom by Vigoureux and 
Grossel @. 

To summarize this discussion at a glance, in Fig. 3 
we have plotted the unit hyperbola [e^ + ^] 2 — [e^-*] 2 = 1, 
assuming real values for all the variables. The interface 
action transforms then the point 1 into the point 0. The 
same hyperbola, when referred to its proper axes, appears 
as g(+)g(-) = 1/2. 



V. THE PHYSICAL MEANING OF INTERFACE 
COMPOSITION 

To conclude, it seems adequate to provide a physical 
picture of the matrix manipulations we have performed 
in this paper. First, the inverse of an interface matrix, as 
pointed out before, describes the interface taken in the 
reverse order. 

Concerning the product of interface matrices, this 
operation has physical meaning only when the second 
medium of the first interface is identical to the first 
medium of the second one. In this case, let us con- 
sider the interfaces 01 and 12. A direct calculation from 
Eqs. (Q) shows that 



loili 



If 



(31) 



for both basic polarizations, which is equivalent to the 
constraints 



'02 



'02 



rpi + rig 
1 + rairu ' 



toitl2 

1 + r Q ir 12 ' 



(32) 



Note that the reflected-amplitude composition behaves 
as a tanh addition law, just as in the famous Einstein 
addition law for collinear velocities: no matter what val- 
ues the reflection amplitudes roi and r\ 2 (subject only 
to | roi| < 1 and |ri 2 | < 1) have, the modulus of the 
composite amplitude |ro2 1 ) cannot exceed the unity. Al- 
ternatively, we have 

r 02 = r 02 = tanh(C 02 /2) - tanh(Coi/2 + & 2 /2) (33) 

which leads directly to the first one of Eqs.(|32|). 

On the contrary, the transmitted amplitudes composes 
as a sech, whose addition law is more involved and is of 
little interest for our purposes here. 

Obviously, for this interface composition to be realis- 
tic one cannot neglect the wave propagation between in- 
terfaces. However, this is not an insuperable drawback. 
Indeed, let us consider a single layer of a transparent ma- 
terial of refractive index N± and thickness d\ sandwiched 
between two semi-infinite media and 2. Let 



2tt 

Pi = —Nidx cos 6^ 
A 



(34) 



denote the phase shift due to the propagation in the layer, 
A being the wavelength in vacuum. A standard calcula- 
tion gives for the reflected and transmitted amplitudes 
by this layer the Airy-like functions H 



-R012 — 



Ti 



012 



rpi + r 12 cxp(-i2fii) 
1 + ?Wi2 exp(-i2/3i) 

t 01 t 12 exp(-z/3i) 
l + ?Wi2 exp(-i2/3i) 



(35) 
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The essential point is that in the limit Pi = 2nir (n = 
0, 1, . . .), which can be reached either when d\ — > 
or when the plate is under resonance conditions, then 
-R012 — * ^02 and T012 — * to2)i and we recover Eqs. (32). 
This gives perfect sense to the matrix operations in this 
work. 



VI. CONCLUSIONS 



in the physics of a plane interface. In these new suitable 
variables the action of any interface appears in a natural 
way as a hyperbolic rotation, which is the natural arena 
of special relativity. 

This formalism does not add any new physical ingre- 
dient to the problem at hand, but allows one to obtain 
previous results (like Fresnel formulas or Stokes relations) 
in a particularly simple and elegant way that appears 
closely related to other fields of physics. 



We have discussed in this paper a simple transforma- 
tion that introduces remarkable simplicity and symmetry 
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FIG. 3: Schematic plot of the hyperbolic rotation performed 
by the interface 01 that transforms on the unit hyperbola the 
point 1 into the point 0. 
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